Exact diagonalization study of the trionic crossover and the trion liquid in the attractive 

three-component Hubbard model 
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We investigate the trion formation and the effective trionic properties in the attractive Hubbard model with 
three fermionic colors using exact diagonalization. The crossover to the trionic regime with colorless compound 
fermions upon increasing strength of the onsite attraction parameter U features smoothly evolving ground state 
properties and exhibits clear similarities to the BCS/B EC-crossover for two colors. In the excitation spectrum, 
there is a clear gap opening between a band of well-defined trions and excitations of broken-up trions at U c ~ 
1.8/. This picture remains the same away from the SU(3)-symmetric point. The spatial pairing correlations for 
colored Cooper pairs are compatible with a power-law at small attractions and change to an exponential decay 
above the trionic crossover. Furthermore, we show that the effective trionic liquid for U > U c can be well 
modeled with spinless 'heavy' fermions interacting with a strong nearest neighbor repulsion. 



I. INTRODUCTION 

The possibility to load several hyperfine species of ultra- 
cold atoms into optical lattices^ has spurred interest in the 
study of multi-color lattice models where the color index is 
used for the internal quantum number of the atoms. In partic- 
ular, fermionic models with three colors were considered^ 
because three colors is the first step beyond the usual two 
internal degrees of freedom of simple spinful many-electron 
systems, and, more ambitiously, these 3-color systems might 
reproduce phenomena known from quantum chromodynam- 
ics where 3 colors of quarks interact and give rise to different 
types of composite particles. Assuming an idealized situation 
where the hopping of the individual colors and the interac- 
tion between them are color-independent, the core model is 
the SU(3)-symmetric Hubbard model with the Hamiltonian 
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with fermionic annihilation (creation) operators c' 1 '. The 
color indexes a, (3 range from 1 to 3. / denotes the hopping 
amplitude between nearest neighbored sites / and j on the lat- 
tice. Note that we have chosen to use the parameter values 
U > for the attractive case. 

In the ultracold atom framework, the strength and sign of 
the onsite interaction U between the atoms can be taken as 
variable. The case of attractive interactions —U < seems 
to offer an appealing trajectory as the strength of the attrac- 
tion is increased. At weak coupling, the two-particle inter- 
action flows to strong coupling at low temperatures, leading 
to the formation of Cooper pairs^^. As the bare attraction 
and the kinetic energy are SU(3) symmetric, there is a de- 
generacy between the energies of all possible two-color pair- 
ing states. In mean field theory for the infinite system where 
symmetry-breaking is possible, this results in a 5-dimensional 
sphere of degenerate BCS ground states, with a vector-like 
order parameter composed out of three complex numbers 
A a p = -A^q., where only the magnitude Aq = YiaB l^ojsl 2 /2 is 



fixed. This state has some rough resemblance with the color- 
superfluid composed out of quark-quark-pairs in high-density 
QCD^. Upon increasing the attraction strength, the formation 
of onsite-molecules of three particles with different colors, so 
called trions, becomes more favorable 4 ' 7 . Indeed, using a vari- 
ational ansatz 4 -, a trionic transition from the color-superfluid 
into a trionic phase with colorless, heavy trions was found. 
The trions are the analogues of hadrons in QCD, and the fa- 
mous QCD transition is supposed to occur between a color- 
ful state at high quark densities and confined hadronic mat- 
ter at lower densities. Interestingly, the QCD phase transition 
is usually described as first order transition, while the trionic 
transition in the attractive Hubbard model was first found to 
be continuous 4 ". Very recently, a refined variational study for 
finite temperatures indicated a first-order transition for the tri- 
onic case as well^,. 

With the gross features of the phase diagram as function 
of U now being clear, many interesting questions arise. Nat- 
urally, one would like to check the results of the variational 
approach 4 with an independent method and learn more about 
the trionic transition. Further the trionic state itself is an in- 
teresting interacting heavy-fermion system and could undergo 
additional phase transitions. However the quasiparticle prop- 
erties and interactions of the trionic state are not easily found. 
One way to proceed, e.g. for finding the dispersion of the 
trions, is a systematic expansion in t/U around the 'atomic' 
limit. Alternatively, the trionic bandwidth and effective inter- 
actions can be addressed using exact diagonalization (ED) of 
small systems which in addition provides a wealth of other 
information, e.g. on the excitation spectrum and correlation 
functions. Here we employ full ED for the calculation of spec- 
tral functions and the Lanczos ED method to study the trion 
formation in the energetically lowest states. Using this infor- 
mation, we derive an effective low energy trion Hamiltonian. 
This Hamiltonian in turn can be analyzed with respect to sub- 
sequent trion phase transitions. 

For the one-dimensional case, additional instabilities of 
three-component systems like density waves or a Fulde- 
Ferrell-Larkin-Ovchinnikov (FFLO) phase have already 
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been discussed, e.g., using density matrix renormalization 
group 7 ' 10 ' 11 . Moreover breached pairing was found for Fermi 
gases in a trap 12 . There also were reports about the successful 
creation of a degenerate Fermi gas consisting of three differ- 
ent hyperfine states of 6 Li (Ref. O as well as mixtures of 



5 Li with 40 K (Ref. Q 



or even 



173 



Yb (Ref. El). Further 



there are theoretical estimates on the stability of such mix- 
tures themselves^ and against particle loss due to three-body 
recombination processes 



U7.18 



II. LOW ENERGY SPECTRUM AND TRIONIC REGIME 

Let us start with the many-particle low energy states of a 
SU(3) symmetric system as found by ED using a Lanczos 
scheme. In order to get a concrete picture, we compute the 
spectrum of a small chain of N sites with periodic boundary 
conditions and n a = 1 fermion per color. In Fig. [T]we show 
the energies of the energetically lowest N + 2 states for weak 
and for stronger attraction. The number of states shown is 
chosen Af + 2, as the interesting trion physics happens in the 
lowest N states, where the (N + l)st state marks the beginning 
of the less trionic regime. For larger interactions, the N lower 
states received their main weight from the basis states with 
the trion on one of the N sites. We find the same sequence of 
energies for N - I trions on N sites, now the trion-hole has 
N possibilities to reside on. In addition, we plot the trionic 
weight w t , defined as the fraction of an eigenstate that is due 
to pure trion states, i.e. those N states with the maximum of 
three fermions on a single site. We observe that, for U = t, 
w, in low energy states is relatively small, and there are strong 
fluctuations in w, for the lowest N + 2 states. Hence the trionic 
weight seems to be a less important property of the low energy 
states. This changes for larger attraction: for U = 6f, there are 
N low energy states which are separated in energy from the 
higher states by a gap A, which grows roughly linearly with 
U minus an offset. For these lowest states, w, is near unity 
and much higher than for the following states. The excited 
states have a lower trionic content and have hence a larger 
amplitude of basis states with broken-up trions. Clearly, this 
data shows the formation of a trionic band, which, for larger 
attraction, is separated from excited states with split-up trions 
by an energy gap of the order \U\. One may also observe from 
the data at large U that the lowest many-particle states have by 
no means the largest trionic weight. Rather the highest states 
in the trionic band are the most trionic. Clearly, this is due to 
a compromise between kinetic and interaction energy, i.e. the 
trions need to split up and sacrifice some binding energy in 
order to lower their kinetic energy. 

With the definition of the trionic annihilation (creation) op- 
erator tf\ via ti = c,ic,2C,3, we take a look at the spatially 
averaged trionic anticommutator |7 T , f] — 1/N \t] ti + tit- ), 
which is a measurement for the well-definedness of the trions 
as compound fermions, and the binding energy of the lowest 
lying N + 2 states. Now the formation of the trionic regime 
becomes even clearer. For weak attraction U - t the lowest 
states in energy have nearly zero interaction energy, that is, the 
expectation value of the interacting part of the Hamiltonian, 
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FIG. 1: (color online). Data for a 14-site-chain, 1 fermion/color. 
Right vertical scale: Empty circles: Energies of the lowest states 
from ED, U = t. Full circles: same for U = 6t. Left vertical scale: 
Trionic weight w, for the same states, empty squares U = t, filled 
squares U = 6r. 



as they basically consist of non-trionic states. We furthermore 
see in Fig. [2] that for U — t the anticommutator does not 
change much over the N + 2 states. It is still close to the non- 
interacting result n 3 + (1 — n) 3 « 0.85 for the density n = 1/19 
per color. In contrast, we find jumps in these values between 
the Mh and the N + 1st state at large attraction U - 6t. Where 
the binding energy of nearly 3U identifies the lowest N states 
to be trionic, the drop to roughly 1 U for the following states 
indicates that these are composed out of a pairs separated spa- 
tially from one single fermion. This suggestion is supported 
by the values of the trionic anticommutator at U = 6t, also 
plotted in Fig. [2] The data clearly shows a drop between 
states N and N + I. 
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FIG. 2: (color online). Data for a 19-site-chain, 1 fermion/color. 
Right vertical scale: Empty circles: interaction energies of the lowest 
states from ED, U = t. Full circles: same for U = 6t. Left vertical 
scale: Trionic anticommutator for these states, empty squares U = t, 
filled squares U = 6t. 
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A closer look at the anticommutator for U — 6t within the 
lowest N states reflects the behavior of the trionic weight de- 
scribed above. While the trion formation is very strong and the 
anticommutator is close to one for all states in the lowest band, 
it is not the groundstate which has the largest anticommutator. 
Instead, from the groundstate to the highest trionic state the 
value of the anticommutator rises monotonically. This can 
again be explained by the fact that the states within the trionic 
band try to optimize their kinetic energy by breaking up the 
trions virtually. 
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FIG. 3: (color online). Data for a 14-site-chain, 1 fermion/color. 
Location of the trion gap opening for asymmetric interactions J7 12 = 
U13 = U > and U23 + U. Inset: Trionic gap vs U/t for the 
symmetric case J/ 12 = I/13 = U23 = U. 

The energy gap between the states with highest trionic 
weight and the rest of the spectrum shall be called the trionic 
gap. In the inset of Fig. [3]we plot the trionic gap 5, as a func- 
tion of U for two different densities. The onset of the gap, and 
hence the trionic band, U c is at 2U C ~bandwidth (4f in this 
case), where 2U is difference of the binding energy of a single 
trion and a 'pair + separated-fermion' configuration. We have 
checked that U c for one trion on N sites does not vary strongly 
with N, and also for more trions it is only changed slightly. 
The value for U c agrees quite well with the location of the 
superfluid-to-trions transitions in the variational treatment of 
Rapp et al., as will be discussed more somewhat below. 

The trionic gap described above shows a marked onset at a 
critical U c . It should be emphasized however that this gap is 
not to be confused with a gap in the single-particle spectrum 
or in the effective low-energy trionic spectrum. It is a prop- 
erty of the excitation spectrum and not of the ground state, 
separating the primarily trionic many-particle states from the 
less trionic ones. Hence the opening of the trionic gap should 
be understood as parameter value where a full band of well- 
defined single trion excitations becomes separated from other 
excitations involving broken-up trions, and where an effective 
trion model which does not resolve the color degree of free- 
dom becomes appropriate. 

Ground state properties evolve completely smoothly as 
function of the interaction strength. In Fig. [4] we plot the tri- 
onic weight in the ground state versus attraction strength U for 



one trion on 10 and 20 sites. The rise of w t contains some steps 
which however are clearly seen to be finite-size-related. As a 
comparison we also plot the weight of doubly-occupied sites 
in the ground states in a system with only two colors and local 
attractions, corresponding to the BCS/BEC-crossover. As the 
local pairs have less binding energy, the crossover occurs at 
larger attraction U, but otherwise the behavior of the weight of 
compound particles in the ground state is very similar. Hence, 
if we do not allow for spontaneous symmetry breaking, there 
is a smooth connection between the weak-coupling ground 
state and the trionic state at large attractions. Nevertheless 
the opening of the trionic gap is a significant phenomenon, 
as it allows us to describe the low-energy theory in terms of 
well-separated trionic excitations only. 
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FIG. 4: (color online). Trionic weight and pairing weight in the 
ground state of 1 1 1 (short for n\ = 1, m = 1, «3 = 1) and 1 10 systems 
for a 10-site-chain and a 20-site-chain. 

In Fig. [3]we also show the location of the trionic gap open- 
ing for interactions which violate the SU(3) symmetry and 
where U23 acting between colors 2 and 3 is different from 
the other two interaction constants. The trion formation is a 
generic feature as long as all three colors attract. A weak on- 
site repulsion between two colors is also tolerable, as one can 
see from the extension of the trionic regime to negative 1/23- 
Hence, while the color-superfluid at weak attraction requires a 
symmetric situation, the trion formation is a general feature of 
the three-color system and could thus be realized much more 
easily in experiments. 

We have also compared the value for the trionic gap open- 
ing for different system sizes with one or two trions and sym- 
metric interactions. We found small quantitative changes, but 
there was always a clear energetic separation between pri- 
marily trionic and other states at interaction strengths above 
U c ~ bandwidth/2. The data suggests clearly that U c remains 
nonzero in the dilute limit. 

An obvious question is whether the system of three colors 
of fermions with attractive interactions also contains a BEC 
regime, where the ground state is composed primarily out of 
strongly bound pairs of two different colors. In the limit of 
local pairs, there is a simple argument to understand that such 
pairs are not stable. If such a difermionic regime would exist, 
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it should be found at [/-values below the trionic regime, i.e. 
at U < U c < 2dt in d dimensions on the hypercubic lattice. 
Now, the binding energy of one trion is 3U, which is 2U more 
than for a pair. So breaking up a trion costs 2U, where the 
kinetic energy gain for the system is about 2dt for the freed 
fermion. Breaking up the pair as well now only costs U, but 
gains Adt of kinetic energy. The latter gain is certainly larger 
if we are below the threshold for trions. Therefore this pair- 
split will be preferred instead of breaking up another trion. 
Accordingly there is no stable difermionic regime and by that 
no BEC regime with such strongly bound pairs. This can also 
be inferred from comparing the curves for the trionic case with 
three colors and the local pair case with only two colors in Fig. 
|4] The local pair regime occurs only deep in the parameter 
range of the trionic regime. 

On the other hand, as discussed in previous works 19 at weak 
U Cooper pairing with loosely bound pairs is possible. The ar- 
gument stated above basically means that the trajectory known 
from two-color system with increasing attraction where the 
pair size shrinks continuously becomes unobservable. Al- 
ready at intermediate U, it is more favorable to form local 
trions, and the non-trionic many-particle states necessary for 
maintaining longer-ranged pairing loose their weights in the 
low-energy sector. 



III. CONNECTION TO THE SUPERFLUID-TO-TRIONS 
TRANSITION AND PAIRING CORRELATIONS 

The variational analysis of Rapp et al4 showed a contin- 
uous quantum phase transition between a color superfluid at 
small attraction to a colorless trionic phase at larger attrac- 
tion at a critical U c ~ 1 .774. This transition is driven by the 
competition between the Cooper pairing instability and the lo- 
cal formation of trions that are implemented as completely 
immobile or frozen in the variational ansatz ofi 4 -. In this tri- 
onic wave-function, the expectation value for any non-local 
pairing correlation (PC) vanishes identically. This forces the 
BCS amplitude to precisely zero when the frozen-trions wave- 
function becomes energetically favorable. Very interestingly, 
in a recent paper Inaba and Sugo 9 * present a variational clus- 
ter approximation study of this transition for the density of 
states of the infinite-dimensional Bethe lattice. They find a 
first order transition between the color superfluid and the tri- 
onic state at low temperatures, hence the precise nature of this 
transition is at least a subtle issue that in theory depends on 
the approximation. 

In our finite systems, there cannot be any spontaneous sym- 
metry breaking. Hence instead of a quantum phase transition 
we rather have a crossover into the trionic regime. Never- 
theless one can expect that the U -values for the formation of 
the well-separated trion band are a good estimate for the at- 
traction necessary to suppress the color superfluid efficiently 
by strongly reducing the weight of non-trionic states in the 
ground state. Indeed, this value, expressed in U over the band- 
width, from the exact diagonalization in one and two dimen- 
sions is very close to the variational results in the limit in infi- 
nite dimensions and to the value in the VCA study 9 . 



The spatial Cooper paring correlations, for simplicity in a 
one-dimensional chain, 

no - x') = ^(c„(jc)c /J (^)cj ; (x')c;(jc'))Gs 



for a + /3 can be readily measured in the ground state found 
by the diagonalization. In the groundstate of a non-interacting 
Fermi liquid one finds a spatial dependence that goes like 



n (A^) 



sin (k F Ax) 
(Ax) 2 



with the distance Ax = x - x'. This gives a node at distance 
Ax = n/kf with the Fermi vector kf . For a system with three 
fermions of each species on a 12-site-chain the node is ex- 
pected at Ax = 4. In the strongly coupled regime one expects 
an exponential decay of the pairing correlations, as doubly oc- 
cupied sites necessary to maintain the pairing are punished by 
a finite energy ~ 2U, and hence the weight of these states in 
the ground state should be suppressed exponentially. The left 
part of Fig. |5]shows the pairing correlations for different U as 
function of the distance. For very small interactions, the node 
of the free solution is clearly visible. In the trionic regime we 
find an exponential decay realized up to Ax — kf/n — 4. For 
larger distances the decay is slower due to finite size effects 
coming in from the other side. The data also shows that for 
intermediate attraction neither of the two descriptions is good 
and one gets a crossover behavior. In fact one would expect 
power-law behavior, but we refrain from determining an ex- 
ponent in our small systems. The main point is the onset of an 
exponential decay above a critical attraction strength. 

The right part of Fig. [5]displays the logarithm of the pairing 
correlations on a distance-vs-interaction grid together with the 
respective contour lines. The solid "vertical" line sketches the 
[/-dependent maximum of the pairing correlations. The data 
shows clearly that a small onsite attraction first enhances the 
Cooper pairing. Above a threshold value below the onset of 
the trion formation, the Cooper pairing is suppressed. This 
behavior of the pairing correlations vs U for our small ID 
system follows closely the variational result for the infinite 
and high-dimensional system of Rapp et al. 4 . Fitting the data 
with a linear combination of the two models 



n(jc) ~ A ■ (e~ fx + e- f{N - x) ) 



+ (l-A)-(U (x) + U (N-x)) 



(2) 



where / is a free fit parameter and < A < 1 interpolates 
between free fermion behavior (A = 0) and exponential decay 
(A = 1). Note that the rise of the weight A of the exponentially 
decaying component is rather steep. It seems well possible 
that in a larger system the rise will become a step. This would 
be consistent with a first order transition. In any case the data 
indicates that the pairing correlations undergo a rapid change 
as function of the attraction strength, making it very plausible 
that the precise nature of the color-symmetry breaking transi- 
tion in the infinite system is a subtle issue. 

We think that the similarity of our data for the small chain 
to the variational results for the thermodynamic limit in Refs. 
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FIG. 5: (color online). Data for a 12-site-chain, 333 filling. Left: 
logarithmic plot of the pairing correlations IT vs. distance d x in the 
groundstate for different U. For very weak interaction ((/ = 0.01, 
open squares) the data shows a dip a k F d x = n as expect from the 
non-interacting case. One can also observe the nearly exponential 
decay for strong interactions {U = 6.00, filled squares). For inter- 
mediate attraction (U = 1.60, open circles) near the crossover to 
trions the pairing correlations are largest. Right: logarithmic map- 
ping of the pairing correlation for distances up to half the system size 
and attractions < U < 5 and contour lines. The solid line shows 
the (/-dependence of the pairing correlations maximum. Top: fitting 
parameter A (with error bars) discriminating the weakly-interacting 
regime from the regime with exponential decay of pairing. Also 
shown are the fit results for the exponential decay factor /, which 
have huge error bars for small interaction where A is near zero and 
the exponential decay plays only a negligible role. 

091 is due to the local nature of the trion formation which is 
to a large degree independent of the dimensionality or lattice 
topology. Nevertheless one might wonder if some kind of 
longer-ranged pairing correlations could persist even in the tri- 
onic regime of the infinite system, as indicated by our fit to the 
pairing correlations that shows a remnant power-law above the 
transition. Hence it is interesting to investigate the superfluid- 
to-trions transition in a framework that goes beyond the treat- 
ment of Ref . 3 by allowing a certain admixture of non-trionic 
states into the ground state that could support Cooper pairing 
even in the trionic regime. The VCA study of Ref. |g| is a first 
step, but it still works in infinite dimensions where the kinetic 
energy of the trions is zero due to the scaling of the single- 
fermion hopping. This could underestimate the admixture of 
non-trionic states to the low energy sector. 

IV. EFFECTIVE TRION HAMILTONIAN 

As the trions are relatively well-defined anti-commuting, 
i.e. fermionic particles for sufficiently large U, we can at- 
tempt to describe the low-energy physics with an effective 
trion Hamiltonian. In principle, the effective trion hopping 
and also the effective interactions can be found analytically 
by a t/U expansion^. The main contributions were already 



given qualitatively in Ref. § Alternatively, with our ED ap- 
proach, the values for hopping parameters can be simply read 
off from the eigenvalues of a small system of N sites with 
one fermion per color, i.e. one trion. The N lowest states for 
U c in the case of a single trion in the system can be associated 
with N momentum eigenstates and nicely fit on a -2f e ff cos k x - 
curve, i.e. nearest neighbor hopping dominates at least well in 
the trionic regime. From this we can read off the bandwidth 
D, = 4f eff of the effective trions band. We find D, = l.5t 3 /U 2 
for sufficiently large U. The third power of t occurs due to 
the three single fermion hoppings necessary to move a trion 
as a whole by one site. As our data are obtained on small sys- 
tems one might worry if a finite-size gap ~ 1 /length spoils 
the detection of the small kinetic energy scale for the trion. 
Note however that the relevant finite-size gaps, measured by 
the difference between neighboring energy eigenvalues within 
the trionic band, also scale down with t 3 /U 2 , i.e. always re- 
main smaller than the effective bandwidth. 

The effects of the interaction between the trions can be 
readily observed in the spatial trion correlation functions 
(n'(x)n'(x')), where n\i) = tjti tells if there is a trion on site i 
or not. As visible in the top part of Fig. |8]the main term is a 
strong nearest neighbor repulsion. The physical origin of this 
is the suppression of a kinetic energy gain of the order t 2 /U 
due to single fermion site fluctuations if two trions sit next to 
each other. This interaction V\ ~ t 2 /U is already compara- 
ble to the effective bandwidth at U c which grows with respect 
to the bandwidth deeper into the trionic regime. Hence the 
trionic liquid is an intrinsically strongly interacting system of 
heavy composite fermions. 

What are the consequences of these interactions? First, near 
half band filling on a bipartite lattice, V\ will favor a density 
wave (DW) ground state of the trions with alternating den- 
sities from site to site. This state is adiabatically connected 
to the density wave state of the colored fermions at weaker 
couplings. In one dimension, this CDW for trions has been 
found numerically by Molina et al^ and Azaria et ali^L. Mov- 
ing the density away from half filling, the DW order will melt, 
possibly with an intermediate regime of phase separation. Yet, 
at one-third filling, it can be seen that V\ still strongly influ- 
ences the trion system even in absence of the DW. The trionic 
Green's function is defined in real space by 

<0|ft|n><n|r,|0> „ <0|*i|m><m|fjl0> 

Gii((jJ) = > h > 

' A ^ to- (E -E„) + 4" o-(E m -Eo) + iif* 

with sums over intermediate states n, m and numerical broad- 
ening r\. The first term corresponds to the emission spectrum, 
i.e. transitions from the ground state into states with one trion 
less, while the second part is the inverse emission spectrum. 
The wavevector-resolved spectral function is the imaginary 
part of the Fourier-transform of G, ; (w) to wavevector k. This 
is shown in Fig. [6] The lower plots in (a) and (b) are for two 
trions in the system. The -2? e ff cos k- trionic band exhibits a 
marked break-up oc t 2 /\U\ between portions concentrated at 
\k\ > n/2 and portions connected to k = above the trionic 
Fermi level. Here the Fermi level is determined by the energy 
where the emission part ends and where the inverse emission 
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part starts, it intersects the band at k — +n/3. The half-filled 
k = +7r/3-states are slightly split up as well. In the upper 
plots of Fig. |6](a) and (b) we show the spectral function for 
just one trion in the system. Now one can see a continuous 
trionic band, but the inverse photoemission part contains an- 
other quasiband shifted upward by an energy oc t 2 /U, which is 
a signature of the interaction of the added trion with the trion 
from the initial state. 




k 
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(b)U/t = 8 

FIG. 6: (color online). Data for a 6-site-chain, 1 (2) fermion/color for 
respective upper graph of a (b). The dashed line marks the chemical 
potential /j. The trionic spectral function shows a cosine band of trion 
states. This band shows in the lower plots in (a) and (b) a gap opening 
above the Fermi level at \k\ = n/2 and some backfolding for the case 
of two trions. This gap (in the lower plot of (a) between u> « -12.5/ 
and a) » -12/) scales with t 2 /\U\ and is larger than typical finite 
size effects of the system measured by energy differences in the band 
between neighbored values. The width of the trion band with the 
gap subtracted scales as t 3 /U 2 . There is an additional split-up of 
the k = ±n/3 excitations that seems to scale in the same way as the 
bandwidth. 

In fact, the many-particle spectrum at low energies with the 
break-up of the trion band can be nicely reproduced with spin- 
less fermions with effective hopping f, and nearest neighbor 
interaction V\ . This comparison allows us to get the value for 



the effective V\ of the trions. We therefore first simulated such 
a spinless fermion system, represented by the Hamiltonian 



H 



{UJ) (U> 



(3) 



with one spinless fermion on 16 sites and compared the width 
of the spectrum AE with the U -dependent bandwidth of the 
according trionic system. 



AE ? 
t eS (U) = — =y— 2 



(4) 



The proportionality factor results to y < 1 .50 where the equal- 
ity holds in the limit of infinitely strong attraction U. With 
this value for the effective trionic hopping we compared the 
gap within the trionic band that occurs due to the strong tri- 
onic interactions for more than one trion in the system with 
the gap in the spectrum of the spinless fermions which be- 
comes A ss V + 2th for large enough V. The second term is a 
kinetic energy loss due to the partial Pauli-blocking of neigh- 
bored sites depending on the density of the spinless fermions, 
h. For a system of two spinless particles on twelve sites this 
gap opens between the 54th and the 55th lowest state. The 3- 
color trionic system develops a gap at the same position, given 
by 



U 

A, = 2— f e ff + 2t eS n,. 



(5) 



By equating the first term to the nearest neighbor interaction 
of the effective trion Hamiltonian, we get the very reasonable 
result Veff = 2yf e ff = 3f 2 /t/ for large enough U. 

In summary, we find, including terms up to t 3 /U 2 , the fol- 
lowing effective Hamiltonian for the trionic sector: 



<U> 



(6) 



(U) 



From these terms we observe that there is a hierarchy of en- 
ergy scales: the dominant term is the nearest neighbor inter- 
action V e ff ~ t 2 /U, where \t/U\ < 1. The trion hopping is 
of order t 3 /U 2 . In this Hamiltonian, higher order corrections 
have not been written. In Fig. Qwe display the convergence of 
trion hopping (or bandwidth respectively) and interactions to- 
ward these values as function of U. Notably, on the 2D square 
lattice the corrections are stronger than in one dimension. 

When we extend our analysis to a 2-dimensional square lat- 
tice, e.g. a square 4x4 cluster, the behavior of the trions re- 
mains qualitatively the same. As the hopping is now allowed 
in more than one direction, the width of the trion band and 
the critical value of U before a trionic gap opens are roughly 
doubled. 

The effective nearest neighbor hopping for the trions ver- 
sus U is shown in Fig. [7] It approaches the same con- 
stant f e ff = l.5t 3 /U 2 for large U. One however observes 
that with increasing U the data in 2D converges from above 
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to this value, while it comes from below in ID. The differ- 
ence is due to an additional higher order hopping contribution 
~ f 5 / U 4 coming from two color fermions hopping directly to 
the neighbored site and 1 color fermion hopping around the 
plaquette. Longer ranged trion hoppings like the one across 
the diagonal also exist, but decay at least as t 6 /U 5 and are 
hence less important. 

effective parameter ID vs 2D 




to- 



5 10 15 

U/t 



FIG. 7: (color online). Data for a 16-site-chain vs 4x4 square cluster, 
111 rilling for effective hopping (left vertical scale), 222 filling to 
find effective repulsion (right vertical scale). Left vertical scale, open 
symbols: the 2D hopping (circles) shows an extra term compared to 
the ID hopping (squares). Right vertical scale, filled symbols: the 
effective repulsion in 2D (circles) indicates the larger U c for the trion 
regime than in ID (squares). 

Fitting the effective hopping dispersion in two dimensions 
with a model including nearest-neighbor (a) , (linear) next- 
nearest-neighbor (b) and diagonal hopping (c), 

t e s(k x ,ky)/t — const 

+ a ■ (cos(fc v ) + cos(^)J 

(7) 

+ b ■ (cos(2£ A .) + cos(2/fc v )) 

+ c ■ (cos(A: v + k y ) + cos(k x - k y )J 

we get the following results for the u — f//f-dependence of 
the parameters 

const — —3m - 3m _1 

a — -3m~ 2 — 12m -4 

, - (8) 

b = -6.7tr 5 - 12 • 6.7m" 7 

c = 2-b. 

Here the prefactors of the second leading order terms in a 
and b are not the result of best fit procedures, but lead to rea- 
sonably small deviation from the numerical data using integer 
numbers and simple fractions. The prefactor 12 in front of the 
£/~ 7 -term in b and especially c - 2b are exact within numer- 
ical precision. This factor of 2 for the diagonal hopping with 
respect to the linear next-to-nearest neighbor hopping can be 



understood easily as there are two ways around the corner to 
the diagonal site. 

As one would have expected, the spatial trion correlations 
on the 2D square lattice exhibit again signatures of the effec- 
tive strong nearest neighbor repulsion. This can be clearly 
seen on the 4x4 plaquette with periodic boundary conditions 
shown in Fig. [8] The nearest neighbor sites (1,0) and (0,1) 
have basically zero probability for both being occupied if 
there is a trion on site (0,0). Considering second nearest 
neighbors, the diagonal neighbors (1,1) and the linear sec- 
ond nearest neighbors (2,0) and (0,2) are equally populated. 
For larger clusters we would expect to find enhanced proba- 
bility for trions on the linear neighbors instead of the diagonal 
ones. This is based on the following idea. As a fluctuating 
trion on site (1,1) has two neighbors (1,0) and (0,1) that are 
affected by Pauli blocking mediated by the fluctuations of the 
trion at (0,0), while on sites (2,0) or (0,2) only fluctuations on 
one neighbored site are partially Pauli blocked. This can be 
confirmed for spinless fermions, but for the trions the largest 
system we can handle here is 4x4, where e.g. (2,0) is only 
one hop away from (-1,0) and is hence also confronted by 
Pauli blocking on two sides, and there is no difference in the 
density correlations on these sites. 

While at first sight, the spatial trion correlations look very 
much like the correlations of spinless fermions with appro- 
priately chosen nearest-neighbor interactions, a detailed com- 
parison shows that for two trions on the small cluster, the sec- 
ond particle is pushed away more strongly than for two spin- 
less fermions. The bottom of Fig. [8] quantifies this differ- 
ence. The expectation value for finding the second fermion 
on site (2,2) in the spinless fermion model is subtracted from 
the normalized value for finding the second trion on site (2,2). 
This difference decreases roughly like I/ -1 , supporting a t-V 
model as effective model for large U, but showing that for 
intermediate U the three-color model is still more complex. 
The best fit of this discrepancy needs an extra U~ 2 term and is 
oc (U~ l - 4i/~ 2 ). Note that longer-ranged effective hoppings 
not included in the spinless-fermion model cannot be the rea- 
son for this difference, as they only start with U~ 5 . 



V. CONCLUSIONS 

With an exact diagonalization study we have identified the 
trionic regime of attractive three-color fermion systems on 
simple one- and two-dimensional lattices. This regime devel- 
ops out of the weak-coupling state with deconfined colored 
fermions through a smooth crossover. The amplitude of non- 
trionic basis states in the ground state becomes smaller with 
increasing attraction, but never reaches zero at any finite U. 
While the ground state properties evolve smoothly, in the ex- 
citation spectrum one finds a gap opening at a critical U < half 
the bandwidth, where a band on trionic single-particle excita- 
tions gets separated from less trionic excitations. This signals 
the onset of a regime with well-defined trions throughout the 
Brillouin zone. We have also shown that the Cooper pairing 
correlations go through a maximum slightly below the trionic 
gap opening when U is increased but get markedly suppressed 
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FIG. 8: (color online). Left Top: Data for a 6-site-chain, 2 fermions/color, U = St. The spatial trionic correlation function shows a strong next 
neighbor repulsion. The value of (;i'(0)n'(0)> < 1 reflects the circumstance that the trions are composite particles that are broken up to some 
degree by quantum fluctuations. Left Bottom: Data for a 4x4 square cluster, 2 fermions/color, U = St. Trionic correlation function shows 
strong next-neighbor repulsion, similar to expectations for spinless fermions. Right: Enhanced expectation value for trions on (2,2) relative 
to expectation for spinless fermions and Independence of enhancement. Best fit uses a (/"'-term plus an additional U '~ 2 -term with negative 
prefactor and has excellent agreement with the data. 



in the trionic regime. This finite-size picture is consistent with 
the previous variational study of the color-superfluid-to-trion 
transition by Rapp et alA On the other hand, our data shows 
that the description of the trionic state should be extended to 
capture the admixture of broken-up-trion states to the ground 
state. These states are necessary for the dynamics of the trions 
and might support remnant Cooper pairing correlations even 
in the trionic regime. In this regime, the trions can be de- 
scribed by an effective model of spinless heavy fermions with 
bandwidth ~ f 3 / U 2 and with a strong nearest neighbor repul- 
sion ~ t 2 /U. We have determined the prefactors for these 
effective parameters. In this sense, the trion liquid is an intrin- 
sically strongly interacting limit wherever the trions are well- 
defined. In the limit of large U, the trions can be described by 
a lattice gas where the particles avoid being nearest neighbors 



and where the kinetic energy is suppressed by a factor of order 
t/U. 

The heavy trion liquid phase itself should be an interest- 
ing laboratory for many-particle physics. As already pointed 
out in the literature, the trions order in a density wave pattern 
on bipartite lattices when there is half a trion per site. Fur- 
thermore, as interacting (quasi-)fermionic particles, the trions 
should undergo superfluid pairing transitions in the odd-parity 
channel at low temperatures, reminiscent of neutron or proton 
superfluid in astrophysics. This will be an interesting topic for 
further research. 

We thank M. Aichhorn, W. Hofstetter, A. Lauchli, J. Ort- 
loff, A. Rapp, C. Toke, and G. Zarand for useful discussion. 
This work was supported by the German Research foundation 
(DFG). 
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